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Externally pressurized gas-lubricated bearings have been used successfully to support 
the shaft of a miniature, high-speed expansion turbine. The bearings incorporate stabilizing 
cavities to suppress whirl at high speeds. The theory of the bearing design is presented in 
detail. 



1. Introduction 

In the field of gas liquefaction and low temperature 
refrigeration, small high-speed expansion turbines 
are often employed as the source of refrigeration. In 
the smaller sizes of machines however, the life and 
speed of conventional bearings are limited, and con- 
sequently a lower limit is set to the practical size of 
such machines. This limitation may be alleviated 
by the use of gas-lubricated bearings to support the 
shaft. The bearings may be lubricated with process 
gas, thus eliminating the risk of contamination with 
oil. Since there is no solid friction, there should be 
no wear and consequently the life of the bearings 
should be of indefinite duration. 

In the operation of a gas-lubricated bearing, the 
surface of the shaft must never touch the surface of 
the bearing, since at the high-surface speeds usually 
employed, such contact may result in seizure and 
the consequent destruction of the bearing, Under 
normal operating conditions the shaft should be 
supported by forces due to pressure in the gas film 
between the shaft and the bearing and these forces 
must be such that the shaft is maintained in a 
position of stable equilibrium. 

In certain types of gas-lubricated bearings, a form 
of instability known as whirl is sometimes observed, 
particularly at high speeds of operation. If the 
speed of the shaft is increased above a critical value, 
the center of the shaft orbits about the center of the 
bearing at a speed approximately one- half the speed 
of the shaft — commonly known as half-speed whirl. 

Whirl may be stable or unstable: Stable whirl is 
characterized by the center of the shaft moving in 
a closed orbit about its former position of equilibrium. 
This is illustrated in figure 1. Unstable whirl (fig. 
2) is characterized by the center of the shaft moving 
along a spiral orbit of increasing amplitude. The 
extent of the increase in amplitude is, of course, 
limited to the extent of the radial clearance between 
the shaft and the bearing. 



The bearing described below has been specially 
designed to ensure stability even at the highest 
speeds of operation with helium gas as the pressur- 
izing medium. The theory of this bearing was pre- 
sented at the First International Symposium on Gas 
Lubricated Bearings [Sixsmith, 1959b]. In the pres- 
ent paper, the theory has been extended and applied 
to the design of the bearings of a helium-expansion 
turbine [Birmingham, Sixsmith, Wilson, 1962]. The 




Figure 1. Whirl: constant amplitude. 




Figure 2. Whirl: increasing amplitude. 
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theoretical approach used in this particular design 
and construction is outlined in the following detailed 
analysis. 

2. Shaft Whirl 

In the diagram shown in figure 2, it is assumed 
that the shaft is vertical so that the only forces 
acting on the shaft are those due to pressure in the 
surrounding gas. At a given instant of time, the 
centers of the bearing and the shaft are situated at 
and 0'. The radial displacement 0' is repre- 
sented by e. The shaft is moving to the left along 
a spiral orbit. The resultant of all the forces due 
to gas pressure is represented by the force vector 
Fx- The instant center of the orbit of 0' is situated 
at 0" . The line 0' 0" bisects the angle between 
Fx and the radius f 0. 

It will be assumed that, over a restricted range, 
the force Fx varies directly as the radius 0', and, 
hence, behaves like a spring. The force Fx may be 
resolved into three components. 

Fxcos 2 0=r(0 / ) 

which pulls the shaft toward the bearing center, 

F x sin e 

which accelerates the shaft center along the orbit, 
and a component 

F x cos B sin 

which provides the energy required to stretch the 
"spring". 

The angular velocity co of 0' about is given by 
the equation 



mu?{0 0') = Y(0 0') = Y* 



or 



-4 



Y 
m 



where m is the effective mass of the shaft in the 
plane of the bearing. 

The linear velocity co(0" 0') of the center of the 
shaft may be resolved into a tangential component 

co(0" 0') cos S = ue 

and a radial component 

co(0" 0') sin 6 = ue tan 6. 

Thus the line of centers may be regarded as rotating 
about the center with an angular velocity co, while 
the radial distance e expands as it rotates. The rate 
of expansion is given by the equation 

e==€o6 a,*tana 

where € is the value of the eccentricity when t=0. 
Thus if the direction of the resultant force lies ahead 




Figure 3. Whirl: decreasing amplitude. 



of the line of centers (0 positive), the diagram ex- 
pands exponentially with time; conversely if the 
direction of the resultant force lies behind the line 
of centers (0 negative) as shown in figure 3, the dia- 
gram shrinks in toward the origin and the bearing 
should be stable, even though the force Fx may be 
quite small. 

In practice, the forces due to the film of gas may 
be far from linear, and at a certain speed of the 
shaft, if the eccentricity is increased from zero, the 
sign of the angle 6 may change from positive to nega- 
tive at some critical value of the eccentricity. Under 
these conditions the diagram should expand asymp- 
totically toward the critical value of the eccentricity 
and whirling of stable amplitude should occur. This 
is often observed in practice. 

3. Fluent Film Lubrication 

The action of a simple sleeve bearing is illustrated 
in figure 4. A shaft of radius r is surrounded by a 
sleeve of radius r+5, the interspace being filled with 




Figure 4. Hydrodynamic pressure distribution. 
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a continuous film of lubricant. The center of the 
shaft 0' is displaced from the center of the sleeve 
by a small distance, e. If the shaft rotates counter- 
clockwise, the film of lubricant is carried into a com 
verging channel on the upper right-hand side of the 
diagram, and out of a diverging channel on the lower 
left-hand side of the diagram. Thus the pressure is 
increased on the upper right-hand side and reduced 
on the lower left-hand side. The resultant force due 
to the pressure distribution is at right angles to the 
eccentricity. This is true, however, only for very 
small values of the eccentricity where compressibil- 
ity effects can be neglected; at larger values of the 
eccentricity, compressibility effects introduce an 
asymmetry and a component of force directed 
toward the center of the bearing is generated. 

Suppose the eccentric displacement of the shaft 
from the center of the bearing is maintained at a 
fixed value by means of a short link free to rotate 
about 0. The forces exerted by the film of gas on 
the shaft will cause the shaft to whirl about the 
center 0. The whirling velocity a; will be the veloc- 
ity at which the gas film is carried into the conver- 
gent channel, i.e., ~ where co, is the angular velocity 

of rotation of the shaft. 
Thus 



01 s 



(i) 



This phenomenon which is known as "half-speed 
whirr' is often encountered in the bearings of high- 
speed spindles, particularly under conditions of light 
loading. 

If the shaft is of mass m, then there will he a 
centrifugal force rnco 2 e due to the whirling which 
acts in the direction of the eccentricity. The 
pressure in the film of lubricant can only exert 
forces at right angles to this direction, and if half- 
speed whirl is occurring, this pressure is zero; 
consequently, if the short link is removed, the shaft 
will move outward along a spiral orbit until it comes 
into contact with the sleeve. As it moves out, 
however, compressibility and other effects cause 
the pressure distribution to become asymmetrical. 
Under these conditions, the new pressure distribution 
provides a force which can be resolved into a tangen- 
tial component and an inwardly directed radial 
component. In practical bearings, small deviations 
from perfection, such as variations in the film 
thickness, may introduce such a radial component 
even when the shaft is at the center of the bearing. 
At moderate speeds this radial component may be 
sufficient to balance the centrifugal force due to 
whirl. 

As the speed is increased, however, the radial 
component of the force due to the film increases 
as the first power of the speed, while the centrifugal 
force due to whirl increases as the square of the speed. 
Thus it is evident that at some critical speed the 
centrifugal force should become equal to the radial 
restoring force, and at this speed the shaft should 



become unstable. This phenomenon is often ob- 
served in practice. 

In a perfectly constructed bearing, however, it 
is evident that a whirl of small amplitude should 
commence even at a very low shaft speed. This 
prediction has been verified experimentally by 
Boeker and Sternlicht [1956]. 

The tangential force F H which is generated in a 
bearing of radius r, of radial clearance 8 and of 
length L, is given by the equation 



F H - 



\%rp.^Lt 



(2+-ft 2 )(l -/r) 



^~°) j '{f") 



(2) 



where 



/x is the viscosity of the gas in lb-sec/in. 2 

e is the eccentricity in inches 

n is the eccentricity ratio 

u s is the rotational velocity 

co is the angular velocity of whirl and 

j I — > n j is the end leakage factor which depends on 
— and n. 



The first, second, and fourth terms on the right- 
hand side of eq (2) are taken from a paper by Ford, 
Harris, and Pantall [1957]. The third term, accord- 
ing to Barwell [1956], is introduced in order to 
account for the effects of whirl. Equation (2) 
holds only where compressibility effects can be ne- 
glected so that the pressure distribution is symmetrical 
about the line of centers. 

The value of the end leakage factor has been de- 
termined experimentally by Ford, Harris, and 
Pantall [1957]. According to the curves in figure 18 
of their paper, the value of the end leakage factor 
in bearings whose length is less than twice the di- 
ameter is approximately 



J \r/ 6.0 



1 L 

r 



(3) 



and for small values of n } the leakage factor is 
independent of the eccentricity ratio. 

Thus, at small values of the eccentricity ratio, the 
tangential force in a sleeve bearing whose length is 
less than twice the diameter is given approximately 
by the equation 



F H =E H e 



2irjj,r 2 L 2 6 /V 



(t-») 



(4) 



where E H is the tangential "spring" modulus. 
Equation (4) may be written 



F H =E H e~ 



2 T 2 

t/jlJ' L co s e 



('-S 



(5) 



The force F H can be either positive or negative 
depending on the value of to. In fact, F H =0 when 
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o) is exactly equal to -—• There is no pressure acting 

on the shaft due to the gas film under this condition. 
The force F H acts at right angles to the displacement 

e, and when co<— > then F H acts in the direction of 

shaft rotation and whirl is accelerated; conversely, 

when o>>-^> then F H acts in the opposite direction 

to shaft rotation and whirl is retarded. 

An appreciation of the magnitude of the tangential 
force in a simple sleeve bearing at an eccentricity e 
may be obtained from the following example. 



Example 1 

Radius of shaft 
Radial clearance 
Length of bearing 
Viscosity of gas 
Angular velocity of shaft 
Angular velocity of whirl 



0.437 in. 

1.50X10" 3 in. 

0.281 in. 

2.97 X10~ 9 lb-sec/in. 2 

1.88X10 4 radians/s 

Zero 



E H - 



3.14X2.97X10- 9 X(0.437) 2 X(0.281) 2 XL88X10 4 



(1.5X10- 3 ) 3 



##=784 lb/in. 



(6) 

(7) 



Thus at an eccentricity of 0.000393 in. (0.001 cm) 
the tangential force is approximately 0.308 lb weight. 

It is evident that in order to insure stability at 
finite speeds of the shaft, the bearing must be de- 
signed so that a radial restoring (centering) force 
F R is generated. If this force is proportional to the 
eccentricity, it behaves like a spring, and it may be 
expressed as 

Ye = F B (8) 

where Y is the modulus of the "spring" and e is the 
eccentricity. An appreciation of the radial spring- 
modulus required in one of the bearings of the 
miniature turbine may be obtained from the follow- 
ing example. 

Example 2 

The effective mass of the shaft in the plane of the 
most heavily loaded bearing is 0.775 lb, and the 
designed rotational speed is 3,000 rps. The corre- 
sponding half -speed whirl velocity is 1,500 rps = 9430 
radians/s. Inserting these values in the equation 
Fe = mco 2 e we get 

F=0.775X9430 2 /32.3X12 (9) 

= 178,200 lb/in. (10) 

At an eccentricity of 0.00039 in., the centrifugal 
force would be approximately 69.5 lb weight. 



4. External Pressurization 

The necessary restoring force may be generated 
by means of external pressurization. The sleeve is 
provided with a number of equally spaced injector 
holes, as shown in figure 5. Compressed gas is 
supplied to the injector holes and flows radially out- 
ward in the clearance space between the shaft and 
the sleeve. The shaft is found to float centrally in 
the sleeve, a well-known effect. 

The reason for this effect becomes evident from an 
examination of a plane circular bearing as shown in 
the insert of figure 6. In this figure the relationship 
between the mean pressure between the plates and 
the plate separation is shown for various supply 
pressures. The graphs in figure 6 show the variation 
of pressure with the separation of the plates for 
plates of % in. diameter. In order to find the 
corresponding variation of pressure for plates of 
any diameter, the scaling relationship h 2 ~D is used 
in conjunction with a constant orifice to pad diameter 
ratio. Accordingly, the graphs shown in figure 7 
are drawn for pads % in. diameter and a plate 
separation of 0.00165 in., which corresponds to a 
plate separation of 0.0015 in. with plates %& in. 
diameter. 

As shown in figure 5, the pads in this bearing are 
not circular but have a rectangular configuration, as 
shown. It will be assumed that each pad is equiva- 
lent to a circular pad % & in. diameter. The mean 
pressure in the clearance space between a pair of 
diametrically opposed pads [Sixsmith, Wilson, Bir- 
mingham, 1961] is shown in figure 7. Curves a and 
b show the variation of mean pressure with clearance 
on the right-hand and left-hand sides of the shaft. 
Curve c shows the variation of differential pressure 
with the eccentric displacement of the shaft. Over 




Figure 5. Conventional gas-bearing arrangement. 
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Figure 6. Load versus plate separation curve. 

a restricted range the differential pressure (curve c) 
may be assumed to be linear, thus the corresponding 
restoring force is also linear, and is directly propor- 
tional to the displacement. 

If the direction of the radial displacement makes 
an angle y with a diameter through the center of a 
pair of diametrically opposed pads, the component 
of the deflection with respect to the pads is e cos 7. 
Similarly, the component of the force generated by 
the pair of diametrically opposed pads in the direc- 
tion of the displacement is also proportional to cos 7. 
Thus the total restoring force along any radius is 
given by 

Fe^^eSLCOS'y]^ 



or 



4 de 



(id 



where Y is the "spring" modulus of the gas film, N 

is the number of pads, A x is the area of a pad, and 

1 dP 

--^— is the differential pressure due to a single pad 

as shown in figure 7. 

Example 3. 

In the present hearing there are eight pads, and 
for this particular number of pads, 



and (11) becomes 



S [cos 2 7^=4 



Ye 



-2A bP e 
-2A X — e. 



(12) 
(13) 




.0010 .0010 

ECCENTRICITY , in. 

Figure 7. Mean pressure on gas-bearing pads. 

The effective area A } of a pad is 0.076 in. 2 , and the 
slope of curve c at 200 psig supply pressure is 44,000 
lb/in. 2 / in. Putting these values in (13) we obtain 
the value of Y. Thus 



F=2X0. 076X44,000 
= 6,700 lb/in. 



(14) 



At an eccentric displacement of 0.00039 in., the 
radial restoring force is 2.61 lb. This is less than 
the corresponding centrifugal force (example 2) by 
a factor of 27. Thus it is evident that the radial 
restoring force, although quite large, is insufficient 
to prevent the shaft from spiralling outward until 
it touches the bearing. It is then necessary to devise 
a method whereby the centrifugal force can be 
reduced to a value below the available radial restor- 
ing force in order to obtain a stable bearing. This 
can be accomplished by reducing the shaft whirl 
speed. The method used to reduce whirl speed is 
presented in sec. 6. 

5. Analysis of Turbine Rotor Instabilty 

In an early design of a turbine [Sixsmith, 1959a] 
in which the turbine was mounted on a flexible shaft, 
it was discovered that a hydrodynamic exciting- 
force can be generated by the turbine rotor. At a 
certain critical pressure of the gas supply, a whirling 
of the shaft would commence. The shaft could be 
made to whirl even when it was prevented from 
rotating. In order to suppress the whirl it was 
necessary to provide damping in the form of oil-filled 
dash-pots. 



720-591—64- 
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h = 8 + e cos 



^F-^ 




Figure 8. Force distribution on turbine rotor. 



An explanation of this phenomenon is as follows. 
The effect appears to be due to forces which are 
generated by the variation of blade-tip leakage with 
rotor eccentricity. Consider the turbine rotor shown 
in figure 8. If the rotor is displaced downwards as 
shown, the tip leakage is reduced around the lower 
half of the rotor and correspondingly increased 
around the upper half. Around the lower half of the 
rotor the tip leakage is less, hence an increased 
thrust is developed on the blades. This force is to 
the right for counterclockwise rotation. Around the 
top half, tip leakage increases and the thrust on the 
blades is reduced. This can also be regarded as an 
added component of force to the right. Conse- 
quently there is a net force to the right acting on the 
rotor as a whole. This force is at right angles to 
the displacement, and therefore, it is in the same 
direction as the hydro dynamic force which is gen- 
erated in a simple sleeve bearing. 

An estimate of its magnitude may be made if we 
assume that the loss of torque due to tip leakage is 
directly proportional to the tip clearance, i.e., 



dT~- 



Td4V 
2tt L 



2wr l hv l 



7rv 2 (r 2 1 —r 2 2 )_ 



(15) 



where 



on the blades 
lb 



-in. 



dT= element of torque acting 

between <j> and <j>-\-d4> 
T= total torque acting on the rotor 
r x =tip radius — inches 
r 2 =root radius — inches 
A=mean tip clearance between <j> and w-\-d<j) — 

inches 
v x = axial component of velocity of gas through 

the tip clearance in./s 
v 2 == axial component of velocity of gas through 

the annulus — in./s between r x and r 2 . 

If the rotor is deflected downwards a distance e 
the tip clearance h at a point which makes an angle 
</> with the direction of the displacement is given by 



h = 8-{-e cos </>. 



(16) 



The horizontal force dF to the right which acts on 
the small segment of blading lying between and 
(j)+d4) is given by 



AW dT * 

dr = — cos 6 

T 



where r m =mean radius of rotor blades. Thus 

dF=—-\ 1 , 2 2 \ cos0. 

27rr m L Tr{rl-rl)v 2 ] 

Substituting (5 + e cos </>) for h we have 

dF= 2^ m \}-J^^ 2 {8+e C ° S 0) J C0S *' 

On integrating between and 2ir, all terms will 
vanish except the one containing cos 2 </>. thus we 
have 



Te 2r 1 v 1 



2wr m {r\—rt)v 2 



G 



sin cos 0+ 



w: 



7rr m (r 2 —r 2 2 )v 2 

Tew 

~r m (r 2 —r 2 2 )v 2 



,(^+i sin2 *) 2 



which may be written as 

F=kTe 



(17) 



where 



Tifh 



r m (r\—r$v 2 



Thus it appears that the whirl-inducing force is 
proportional to the product of the torque and the 
eccentricity, but is independent of tip clearance as 
long as the assumption that the blade passage 
efficiency is inversely proportional to tip clearance is 
valid. Since for a given eccentricity it is propor- 
tional to the torque, it is capable of inducing a 
whirl even when the shaft is prevented from rotating. 

In the miniature turbine the value of the torque 
at 3000 rps and 3 kw of power is 1.406 lb-in. The 
values of r lt r m , and r 2 are 0.436, 0.397, and 0.36 in. 
Since the blade outlet angle is 30°, the value of 



we have 



is 2, approximately. Putting these values in (17) 



1.406X0.436X2X6 
z 0.397(0.436 2 -0.360 2 ) 



or 



F=51.15elb. 



(18) 



Thus at an eccentricity of 0.00039 in. the tangential 
force is 0.02 lb weight. 



106 



RESISTANCES R 



CAVITIES OF 
CAPACITY C 




/^■mr *\ 



Figure 9. Gas hearing with stabilizing cavities. 



6. The Analysis of a Method of Stabilizing 
a Gas-Lubricated Bearing 

In order to ensure stability, it is necessary to reduce 
the angular velocity of whirl until the centrifugal 

force becomes less than the radial restoring force. 
A reduction in the angular velocity of whirl may be 
effected if the bearing is designed to generate a force 
which acts in the opposite direction to the hydro- 
dynamic exciting forces due to the bearing and the 
turbine rotor. 

A bearing designed to produce this stabilizing 
force is shown in figure 9. In addition to the injector 
holes, as shown in figure 5, the shaft is surrounded 
by a number of pneumatic phase shift networks 
which are spaced midway between the injector holes. 
Each network consists of a cavity which acts as a 
pneumatic capacitor and a leak orifice which acts 
as a pneumatic resistor. A stream of gas flows into 
each cavity from the adjacent injector holes, and 
escapes from the cavity via the leak orifice. The 
pressure in the cavity is determined by the flow rate 
and the resistance of the orifice. In the event of a 
whirling of the shaft, the flow of gas into each cavity 
is modulated, the modulation being in phase 4 with 
the eccentricity vector. As a consequence of the 
modulated inflow, the pressure in the cavity is like- 
wise modulated, but by virtue of the finite volume 
of the cavity, the phase of the pressure is lagging 
with respect to the phase of the current and the 
eccentricity. The force on the area of shaft facing 



the cavity, being in phase with the pressure, is like- 
wise lagging with respect to the current and eccen- 
tricity vectors. This force can be resolved into 
components whose directions are radial and tan- 
gential with respect to the eccentricity vector. 
The tangential component acts in the opposite 
direction to the hydrodynamic force due to the 
bearing and the turbine rotor. By suitable design 
of the bearing, the tangential component of force 
can be made to exceed the exciting forces so that 
the net tangential force is opposed to the tangential 
motion of whirl, and consequently, the whirling 
motion should be damped out. 

The radial component acts in the opposite direc- 
tion to the centering force due to the injector holes, 
and therefore it reduces the load carrying capacity 
of the bearing. This is inherent in the design and 
cannot be avoided. 

As shown in figure 9, a shaft of radius r is sur- 
rounded by a sleeve of internal radius r-\-8. Gas from 
eight injector holes each of radius r,, equally dis- 
placed around the circumference of the sleeve, flows 
into eight cavities each of capacity C. The ma- 
jority of the gas escapes from each cavity through a 
leak orifice of resistance R. The remainder escapes 
at the ends of the bearing, but it will be assumed that 
this portion is negligible. Each cavity and leak 




Figures 10-15 (inc.) Vector diagrams illustrating the action 
of a stabilizing cavity. 
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behaves like an electrical resistance-capacity net- 
work : the capacity of a cavity being defined by the 
equation 



n bP jr dp T7 p dP T 1U/ 



and the resistance of a leak by 
~=Il lb/s 



(19) 



(20) 



where V is the volume of the cavity, p is the density of 
the gas in the cavity, P is the pressure in the cavity, 
I c is the current of gas accumulating in the cavity, 
and I L is the current flowing through the leak: t is 
time. 

A cavity receives its gas supply from both of the 
adjacent injector holes. Suppose now that the 
shaft is whirling with an eccentricity e and a counter- 
clockwise angular velocity co. The motion of the 
center of the shaft may be represented by a complex 
vector e rotating counterclockwise with the same 
angular velocity (fig. 10). 

Referring again to figure 9, it will be seen that 
the current of gas 1 entering cavity 1 will be given by 



7=iy2/ 1 +l/2/ a 

=Tr j vp l (h l +h 2 ) 



(21) 



where v is the velocity of the gas at the edges of the 
injector holes and p T is its density. The values of 
the clearance at the injector holes adjacent to cavity 
1 are given by 



and 



A 1= =5-f e cos ( <^~Tf) 

h 2 = 8 + e cos ( cot-\- jf J 



where N is the number of cavities. 

The time t is assumed to be zero when e is pointing 
vertically downward through the center of cavity 1. 

Equation (21) may now be written 

I=rr j vpA 8 + e cos (cot — ~J+8 + e cos f <*>t+j^J > 



=irr j vpi*i 2<5+2e cos cot cos 



N 



}■ 



(22) 



The current flowing into cavity 1 thus consists of 
a steady component i which is given by 



I = 2Trr j vp l 8 



(23) 



and an alternating component Ij which rotates in 
phase with the eccentricity vector e and whose 



instantaneous value is given by 

Ij cos o)t=27rr j vpie cos cot cos 



N 



(24) 



In the above analysis, the effects of shaft rotation 
and whirl have been neglected. These effects influ- 
ence the flow rate into the cavity. Referring to 
figure 9, it is seen that the clearance between the 
shaft and the bearing on the right-hand side of the 
cavity is greater than the clearance on the left-hand 
side. As a consequence of the rotation of the shaft, 
more gas is carried out of the cavity on the right- 
hand side than is carried in on the left-hand side. 
Thus there is a net flow of gas out of the cavity. An 
estimate of the value of this component of current 
may be made if it is assumed that the principle of 
superposition may be applied; i.e., that the velocity 
of the gas at any point in the film is the vector sum 
of the radial velocity from the injector holes and the 
tangential velocity due to hydrodynamic action. 

The instantaneous value of the hydrodynamic 
component of current I a is given by 



& s 



7 a sin cot=~L c rp 2 (h2—h 1 ) 



=^L c rp 2 4 8+ecos(o)t+~ 
-[5+6 cos (^-^)]} 

= — co s L c rp 2 € sin cot sin 



N 



(25) 



where l a is the alternating flow of current, o s is the 
angular speed of the shaft, ris the radius of the shaft, 
L c is the effective length of the cavity, and p 2 is the 
density of the gas in the clearance space. 

Referring again to figure 9, it is seen that as a 
consequence of the whirling motion, the shaft is 
approaching the cavity and therefore the volume of 
the cavity is decreasing. This may be regarded as a 
current I w flowing into the cavity, the volume of 
the cavity being regarded as constant. The in- 
stantaneous value of the current is given by the 
equation 



I w sin cot=2coL c rp 2 6 sin cot sin -, 



N 



(26) 



where I w is the alternating flow of current and co is 
the angular velocity of whirl. 

Thus the resultant current I H due to the effects ol 
shaft rotation and whirl is given by the equation 



r T I 2co\ 

I H sin cot =—L c rp 2 eco s ( 1 J 



sm co 



sin^ (27) 



This current may be represented by the vector Ih 
(fig. 11) which rotates 90° ahead of the eccentricity 
vector e. 
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The resultant current flowing into the cavity is 
represented by the vector F. It is seen that 



1' = /, sec 0. 



(28) 



The direct current I flowing into the cavity may 
he represented by vertical line OA, figure 11, and 
the instantaneous value of the alternating current 
by the vertical line OB. This is the projection on a 
diameter through cavity 1 of the alternating current 
V which rotates ahead of the eccentricity vector e 
with a phase angle /3. 

The value of J8 is given by 



tan p=-f-= 



I l{ L c rp 2 a> s 



2wr j vp l 



0-S) 



tun 



N 



(29) 



The steady current I flowing through the resist- 
ance R generates a pressure given by 



/v-io/.'. 



(30) 



This is represented by the vertical line 01), figure 12. 
The alternating current generates an alternating 
pressure component 



P'=rZ=IjZ sec 0= 



TjR sec g IjR sec fl 



s/l + u 2 C 2 R 2 Vl+tan 2 (<*+|(3) 

(31) 



where Z is'the impedance of the resistance-capacity 
network. With respect to the alternating component 
of current, the capacity and resistance are effectively 
in parallel and, consequently, the phase of the 
pressure lags behind that of the current and the 
eccentricity. The phase angle a is given by 



tan (a + p)=wRC. 



(32) 



The currents due to the alternating pressure P' 
are represented vectoriallv in figure 13. The vector 
Ir represents the alternating current flowing in the 
resistance B, and the vector Ic represents the alter- 
nating current flowing in the capacity. The values 
of I R and I c are given by the equations 



1 



?1 
R ~r 



and 



l = p' 03 C=P'coV 



I' 



(33) 
(34) 



The resultant of these two currents must be equal 
to and in phase with the alternating current I' 
flowing into the cavity. 

The pressures P and P' acting on the effective 1 
area A 2 of the shaft, which forms one wall of the 
cavity, generate forces 



F =A,P 



(35) 



and 



F'=A 2 P'. 



(36) 



These forces may be represented by the vertical 
vector OG (fig. 14) and the alternating vector OJ, 
respectively. The alternating vector F' may be 
resolved into a component F' cos a which is in 
phase with the eccentricity and a component F' sin 
a which is lagging 90° behind the eccentricity (fig. 15). 

These forces are generated by a single cavity and 
act along the radius from the center of the shaft 
through the center of that cavity. With a sym- 
metrical array of cavities, the resultant of the steady 
forces Fq will be zero ; this, however, is not true of the 
alternating components. 

The resultant force on the shaft at any instant is 
found by adding vcctorially the instantaneous values 
of the forces due to each of the cavities. 

The instantaneous value of the force due to cavity 
1 is 

F' cos (at-a). (37) 

This is the projection of the alternating vector F' 
along a diameter through the center of the cavity. 
At a given instant it is represented in magnitude and 
direction by the vertical line ON as shown in figure 
14. 

The polar diagram of the force due to each of the 
eight cavities is shown in figure 16. At an instant 
of time t, the vectors numbered 1 to 8 represent the 
forces due to the corresponding cavities. 

The polar diagram represents a field of force ro- 
tating with the angular velocity of whirl. This is 
illustrated in figure 17, which shows the values of 
the eight forces due to the cavities when t has in- 
creased to t+At. These forces can be combined into 
a resultant which by symmetry will lie along OJ. 
The magnitude of this force may be written as F s 




Figures 16 and 17. Polar diagrams of force generated by 
stabilizing cavities. 



i Experimental testing has shown A/2 to be a conservative estimate for "the 
value of A2 as illustrated in figure 20. 
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and is given by 



N 



(38) 



Substituting from (36), (31), (29), and (24) this 
becomes 



or 



where 



F s =NirA 2 r j vpiZe cos -^ sec fi 



F s =GZe 



G=NirA 2 r j vpi cos -^ sec /3. 



(39) 
(40) 
(41) 



G is a constant which depends on the design of the 
bearing and the speed of the shaft. The alternating 
vector GZe may be resolved into components GZe 
cos a in phase with the eccentricity vector or GZe 
sin a lagging 90° behind the eccentricity vector. 
The component of force GZe cos a due to the stabi- 
lizing cavities is in the same direction as the centrif- 
ugal force due to whirl and in the opposite direction 
to the centering force due to the injector holes 
defined in (11). Consequently, the centering force 
is reduced to Ye— GZe cos a. Whirl can set in at 
an angular velocity given by 



i.e., 



mo) 2 €= Ye— GZe cos a; 
/ Y-GZ COS a \* 



(42) 
(43) 



These forces are represented as vectors in figure 18. 




The component GZe sin a provides a force at 
right angles to the eccentricity vector and acts in 
the opposite direction to the forces due to hydro- 
dynamic action or excitation in the bearing and the 
turbine rotor F (17). This is the damping force 
which is required to provide stability. Provided it 
is greater than the sum of the bearing and turbine 
rotor excitations, the direction of the resultant of 
all the forces lies behind the line of centers, and any 
whirl, such as might be excited by an external 
shock, is attenuated until its amplitude becomes 
zero. 

An estimate of the magnitude of the hydrodynamic 
force due to the bearing may be made as follows. 
The bearing is regarded as a simple sleeve bearing 
in which the shaft rotates with an angular velocity 
a> s and a whirl velocity co. The value of the whirl 
velocity co is given by (43). The hydrodynamic 
force which is generated at an eccentricity e can be 
calculated with the aid of (5). The force thus 
calculated may be expressed in the form 



E H e=F H 



(44) 



for small values of the eccentricity ratio e/8. E H is 
a constant whose value depends on the design of the 
bearing. In the above calculation no account is 
taken of the disturbance in the flow pattern due to 
the various holes and slots in the surface of the 
bearing. These have the effect of reducing the 
pressures due to hydrodynamic action and therefore 
the calculated value of the hydrodynamic force E H e 
is an overestimate. Thus it errs on the safe side. 
The force E H e is represented vectorially in figure 18. 
The rotor excitation kT is given by (17). This 
applies however to the force in plane of the rotor; 
hence, the value must be multiplied by some factor 
k\, depending on the overhang, in order to give the 
force in the plane of the bearing. This force may 
thus be expressed in the form 



where 



F T =k'Te 



(45) 
(46) 



This force is shown vectorially in figure 18. The 
term k'T can be regarded as the rotor excitation 
reduced to the plane of the bearing. Keferring to 
figure 19, it is seen that the excitation clue to the 
rotor is greater in the inner bearing than it is in 
the outer bearing. For the inner bearing the value 
of k\ is given by 

Kefs) (47) 

The stabilizing force GZe sin a which is generated 
by the bearing must be greater than the sum of the 
two excitations outlined above if the shaft is to be 
stable. 

The condition required for dynamic stability is 
thus given by 



Figure IS. Illustration of forces on shaft. 



GZe sin a^>k' Te + Ene. 



(48) 
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Figure 19. Shaft resonant system. 

Provided this condition is satisfied, the magnitude of 
e will progressively decrease to zero and the vector 
diagram in figure 18 will shrink into the origin. 
On the other hand, if the sum of the excitation is 
greater than the damping term, the vector diagram 
will progressively expand until e becomes equal to 
the radial clearance — at which point the bearing 
may be destroyed. 

In the above 4 analysis the conditions which deter- 
mine the dynamic stability of the bearing have been 
outlined. However, the ben ring must also possess 
static stability, i.e., it must be capable of supporting 
a radial load, such as the weight of the shaft if it 
were horizontally supported. In this case, co becomes 
zero and if co, is also zero then is zero, and since tan 
(a J rl3)=ojGI\, the phase angle a also takes the value 
of zero. The stabilizing force GZe sin a thus reduces 
to zero, and the decent ering force GZe cos a reduces 
to the value GRe. The centering force now becomes 



F A =(Y-GR)e, 



(49) 



and it is evident that GR should be less than )' in 
order to ensure static stability. If GR should be 
greater than Y, the shaft will be pushed to one side 
of the bearing. It is evident that the value of the 
leak resistance must be carefully chosen; it must be 
large enough to ensure dynamic stability; if it is 
made too large, however, the shaft will not float. 
If we write 

GR 



=B, 



(50) 



then the static stability condition will be satisfied if 
B is less than one. The static centering force may 
now be written 



In the present bearing, B was chosen to be 0.75, 
and consequently, the static load carrying capacity 
is 0.25 Ye. The area of the shaft A 2 which faces each 
cavity is estimated to be as shown in figure 20. The 
value of /3 is given by (29). As a first approximation 
it may be assumed that (3 = and sec /3=1. The 
(low velocity v is assumed to be sonic at the throat 
formed by the annulus 2icrfi\ i.e., the discharge 
coefficient is assumed to be unity. The diameter of 
the leak orifices may be calculated with the aid of 
this equation. The discharge coefficient for a 
circular orifice is approximately 0.6. 

Since the shaft is supported on two bearings, each 
of which behaves Like a spring, it has two principal 
modes and frequencies of whirl. The calculation of 
the modes and frequencies of such a resonant system 
is described in standard textbooks on mechanical 
vibrations. However, for the convenience of the 
reader, the procedure in this specific application is 
outlined below. 

The shaft of the miniature turbine and its supports 
are shown schematically in figure 19, where the 
following symbols and the corresponding values are 
used. 

m=Mass of shaft — lb mass 

r ff =Radius of gyration of shaft about the center of 

gravity — in. 
Yi= Stiffness of gas film in outer bearing— lb/in. 
F 2 z=Stiffness of gas film in inner bearing — lb/in. 
Li=Distance of outer bearing from center of 

gravity in. 
L 2 = Distance of inner bearing from center of 

gravity — in. 

It is assumed that the two radial spring constants 
arc equal, i.e., 



Y l = Y 2 = Y-GZ cos a. 



(54) 



In this equation, a value for a which should 
provide the optimum damping is required. This 
value is given by the equation 



90° -0 



where 



tan 0= 



Ih 



(55) 
(56) 



F A = Ye(l-B). (51) 

The value of the leak resistance is given by 

(52) 



R= P ' = BY < 



Ir Gt 

Substituting for G from (41) we get 

BY 



E= 



NirA 2 rjVPi cos t= sec /3 



(53) 




Figure 20. Effective area facing cavity. 



in 



as shown in figure 11. But I H cannot be evaluated 
until the value of co is known, and this depends on 
the value of cos <x. 

This difficulty may be circumvented if it is as- 
sumed as a first approximation that the angle fi is 
zero and the angle a is 45°. With these values of 
a and fi, Z becomes equal to R/^J2 and (54) becomes 



Y l =Y 2 =Y-}- GR. 



(57) 



In the present calculation the value of R is 
arbitrarily chosen so that 



GR=$ Y. 

4 



(58) 



The two frequencies of oscillation are given by 
the equations 



1 



(a+c) + ^\(a-cy+(^ 



(59) 



and 



where 



2 
00 2 = 



\ («+<0-Vi («-<0 2 +(*) 2 (60) 



7 L2I 2 -^1-*1 

= ; 

m 

■^ll 1 I L 2 1 2 



c= 



mrt 



Equations (59) and (60) are approximate and do 
not include the gyroscopic effects associated with 
rotation of the shaft. However, such effects can be 
ignored for the case considered here. 

The distances of the nodes from the center of 
gravity are given by the equations 



u 



coi — a 



L 4 = 



u%—a 



(61) 



(62) 



The positions of the two nodes are illustrated in 
figure 19. From the positions of the nodes it is 
seen that the predominant amplitudes of whirl in 
the inner and outer bearings occur at the lower and 
higher frequencies of whirl. 

These values of o^ and co 2 may now be employed 
in (31) and (54) in order to find more precise values 
of Y t and Y 2 . From the revised values of Y x and 
Y 2 , a more refined calculation of the values of o? x 
and co 2 may be made. However, in view of the many 



uncertainties in the computation, the first approxi- 
mation is usually sufficiently accurate for practical 
purposes. 

The value of the phase angle /3 may now be 
calculated with the aid of (29). The optimum value 
of the phase angle a is given by (55) . 

The values of the stabilizing cavities may now be 
calculated. We have 



tan (a+P)=a:CR= 



oVpR 



or 



and 



Vi= 



tan (a+/3)P 



V 2 



COipR 

tan (a+P)P 



oo 2 pR 



(63) 
(64) 
(65) 



where cox and co 2 are the whirl angular velocities of 
the outer and inner bearings, V x and V 2 are the 
corresponding volumes of the stabilizing cavities, and 
p is the density of the gas at P units of pressure. 

Referring to figures 12 and 11, it is seen that the 
vectors of pressure and current are given by the 
equations 

P'=rZ (66) 

and 

r=Ij sec p. (67) 



The value of Z is given by the equation 

J? 

Z= {l+tan 2 (a+/3)} w " 



(68) 



The value of the force F' , figure 14, due to a 
single cavity is given by the equation 



F'=AoP' 



(69) 



where A 2 is the effective area of the shaft (fig. 20) 
which forms part of the walls of the cavity. The 
value of the force GZe (fig. 18) is given by the 
equation 

GZ*~F'. (70) 



Thus the value of GZe is given by 



GZe 



NA 2 RIj sec /3 

= 2[l+tan 2 (a:+^)} 1/ 



(71) 



The corresponding tangential and radial com- 
ponents of force are GZe sin a and GZe cos a. 

The condition for dynamic stability is given by 
(48); i.e., 

GZe sin a>k'Te+E H e. (72) 

Provided this condition is satisfied, the direction of 
the resultant of all the forces acting on the shaft 
lies behind the line of centers as shown in figure 18. 
Thus the orbit of the center of the shaft should be 
a spiral of decreasing amplitude as shown in figure 3. 
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l - BY6 sec/3(i- sin/3) 




The values of Y x and Y 2 according to (57) are 
given more simply by the equation 



Figure 21. Vector diagram of forces illustrating design 
procedure. 



Consequently, the bearing should be stable. This 
was proved experimentally in the bearings of the 
miniature turbine which is described by Birmingham, 
Sixsmith, Wilson [1962] 

The design procedure and the relationship between 
the various forces may be illustrated with the aid 
of the geometrical construction shown in figure 21. 
The angle fi is preferably chosen to be less than 
20°, assuming zero whirl velocity, and the diameter 
of the injector holes is calculated according to (29) 
to satisfy this condition. 

In figure 21 we have: 



0'S=BYe 

0'D=BYesec(3 

AN=l/2BYesec 

0'H=AN-0' A sin (3 



(73) 
(74) 
(75) 



= l/2BYe sec /3(l-sin /3). (76) 

Thus the stabilizing force O'H may be expressed 
directly in terms of the centering force Ye. The 
condition for dynamic stability may now be expressed 
as 

\j2BYe sec 0(1 -sin P)>k'Te+E H e. (77) 

In a trial design the number of pads and the 
radial clearance are adjusted until the above sta- 
bility condition is satisfied. Then the diameter of 
the leak orifices and the volume of the stabilizing 
cavities may be calculated as follows. 



¥,=¥,=¥ 



H*> 



(78) 



and then the values of the whirl angular velocities 
may be calculated with the aid of (59) and (60). 
The diameter of the leak orifices may be calculated 
according to (53) . Then the values of the stabilizing 
cavities may be calculated according to (64) and 
(65). This completes the design. 

In the theory outlined above, it is assumed that 
the leakage along the shaft from the stabilizing 
cavities is negligible. Such leakage introduces a 
phase shift which tends to promote instability. 
Conversely, a leakage along the shaft and into the 
stabilizing cavities tends to promote stability. 
Consequently, it is desirable that the pressure at 
the ends of the bearing should be equal to or greater 
than the static pressure in the stabilizing cavities. 



7. Conclusion 

It is possible and quite practical to provide stable, 
externally pressurized gas-lubricated bearings for 
high-speed shaft operation. The stabilizing cavities 
employed in this design are a definite asset to a 
bearing system of this type and have been shown 
by actual tests to suppress shaft whirl inherent in 
externally pressurized bearings. 

The bearing described here is adaptable to small 
shafts (one inch diameter and smaller) which are 
lightly loaded and require speeds well beyond the 
allowable operating range of rolling contact bearings. 
For miniature applications, the size of the bearing 
and shaft system is limited only by practical fabri- 
cation techniques. 

This type of bearing is particularly advantageous 
for miniature turbine expanders used in cryogenic 
refrigerator systems, where the process gas can be 
used as the bearing gas medium. In this arrange- 
ment, the problem of contaminating the process gas 
from the bearing lubricant is eliminated. It also 
eliminates a mechanical seal between the bearings 
and the process expansion system. 

A prototype turbine expander, incorporating ex- 
ternally pressurized gas-lubricated bearings, was 
fabricated and tested with helium gas as the operat- 
ing medium. The bearings were designed using the 
criterion described in this manuscript. The tests 
showed that the bearings were reliable 4 and performed 
quite satisfactorily. No detectable shaft whirl was 
encountered up to the design speed of the turbine 
shaft. 

8. Summary of Nomenclature 



^4.!= Area of pad — in. 2 

A 2 = Effective area of shaft surface 

cavity — in. 2 
B = Constant defined in eq (50) 



facing 



720-591—64- 
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Ail 



= Capacity of cavity 
= Diameter of pad — in. 
= Tangential "spring" modulus — lb/in. 
= Horizontal force on blading — lb 
= Static centering force — lb 
= Force due to gas pressure on shaft — lb 
= Tangential force — lb 
= Radial restoring force — lb 
= Steady force — lb 
= Alternating vector of force — lb 
= Vector sum of forces due to N cavities — lb 
= Constant defined in eq (41) 
=Meantip clearance between 4> and 0+rf</> — 
in. 

[ = Clearances between shaft and bearing. 
See figure 9 — in. 

= Alternating component of current flowing 
in cavity — lb/s 

= Current through leak — lb/s 

=Flow of gas into cavity — lb/s 

= Steady component of flow into cavity— 
lb/s 

= Alternating component of flow into cav- 
ity — lb/s 

= Alternating flow due to shaft rotation (eq 
25) 

= Alternating flow due to shaft whirl 

= Alternating flow due to the combined effects 
of shaft rotation and whirl 

= Alternating flow into cavity 

^Constant defined in eq (17) 

= Correction for overhang of rotor 

= Constant defined in eq (46) 

= Length of bearing — in. 

= Effective length of stabilizing cavity — in. 



= See figure 19 



stabilizing cavities 

3 



in outer 



^=Mass of shaft in plane of bearing- 
T= Number of pads in bearing 

i=-= Eccentricity ratio 



-lb 



= Center of bearing 
0' = Center of shaft (see fig. 1) 
0"=lnstant center of orbit of shaft (see fig. 2) 
P = Pressure in cavity — psi 
P = Steady component of pressure in cavity — 

psi 
P f = Alternating component of pressure in 

cavity 
R= Resistance of leak orifice — sec/in. 2 
r= Radius of shaft — in. 
Vi = Turbine rotor tip radius — in. 
r 2 = Turbine rotor root radius — in. 
r TO =Turbine rotor mean radius — in. 
r,=Radius of gyration of shaft about its 

center of gravity — in. 
T= Total torque acting on turbine rotor — in. -lb 
t=time — sec 



V= Volume of stabilizing cavity (see eq (19)) — 

in. 3 
Vi= Volume of 

bearing — in: 
V 2 = Volume of stabilizing cavities in inner 
bearing — in. 3 
v= Velocity of gas at edge of injector hole — fps 
Vi= Axial component of velocity of gas through 

rotor tip clearance — fps 
v 2 = Axial component of velocity of gas through 

the annulus between r Y and r 2 — fps 
Y= Radial spring modulus — lb/in. 
Yi = " Spring" modulus of gas film in outer 

bearing — lb/in. 
Y 2 = "Spring" modulus of gas film in inner 

bearing — lb/in. 
Z— Impedance of resistance capacitance net- 
work 
Greek 

a = Phase angle between complex pressure 

vector and eccentricity e 
P= Phase angle between resultant current I r 

and eccentricity e 
7= Angle between e and a diameter through the 
center of a pair of diametrically opposed 
pads 
8 = Radial clearance — in. 
e= Shaft eccentricity — in. 
0=See figure 2 

fi= Viscosity of gas — lb-sec/in. 2 
p = Density of gas in cavity — lb/ft 3 
P!= Density of gas at edge of injector hole — 

lb/ft 3 
p 2 = Density of gas in clearance space between 

shaft and bearing — lb/ft 3 
4>=See figure 8 

co = Angular velocity of whirl — radians/s 
a?! = Angular velocity of whirl in outer bearing — 

radians/s 
co 2 = Angular velocity of whirl in inner bearing — 

radians/s 
co s = Angular speed of shaft — radians/s 
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